Molecular Geometry, Vibrational Assignments, HOMO-LUMO, Mulliken’s charge analysis and DFT Calculations of 2-(2-Phenylaminothiazole-5-oyl)1- methyl-6-methylbenzimidazole

 

S. Sangeetha1*, T. F. Abbs Fen Reji2

1Department of Chemistry, Sivanthi Adithanar College, Pillayarpuram-629501, Tamilnadu, India

2Department of Chemistry and Research Centre, Nesamony Memorial Christian College, Marthandam-629165, Tamilnadu, India

*Corresponding Author E-mail: sangeethasss1982@gmail.com

 

ABSTRACT:

The vibrational wave numbers of 2-(2-phenylaminothiazol-5-oyl)-1-methyl-6-methylbenzimidazole were calculated using B3LYP/6-31G basis set and calculations are used to assign vibrational bands obtained experimentally. The B3LYP method is able to predict vibrational frequencies and structural parameters. The optimized molecular geometry, bond lengths, bond angles, dihedral angles and harmonic vibrational wave numbers of the titled Compound have been investigated by Density Functional Theory (DFT) method. The geometries obtained from DFT method is found to be in good agreement with experimental data. The Mulliken population analysis on atomic charges has been computed using DFT calculations. Energetics of the Highest Occupied Molecular Orbital (HOMO) and Lowest unoccupied Molecular (LUMO) of the molecule were calculated using the Gaussian 09 software package.

 

KEYWORDS: Benzimidazole, B3LYP, Vibrational frequency, DFT, HOMO, LUMO, Gaussian.

 

 


INTRODUCTION:

Heterocyclic compounds have played an important role in pharmaceutical chemistry due to their biological activities1. Benzimidazole shows several structural features and biological activities just an indoles2. The literature survey shows several examples of compounds having benzimidazole derivatives possess anticancer3, anti-inflammatory4, antibacterial5 and anti-fungal6 activities.

 

Vibrational spectroscopy is used extensively in organic chemistry for the identification of functional groups of organic compounds, for studies on molecular confirmation, reaction kinetics etc7-9.

 

 

 

Literature survey reveals that to the best of our knowledge, the results based on quantum chemical calculations, HOMO-LUMO analysis of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole have no reports. Here in we reported detailed interpretations of vibrational assignments which are acceptable and supportable.

 

Patil et al. reported the DFT study on dihydroxyphenylbenzothiazole by using B3LYP/6-31G(d)10. The main objective of this paper is to present, more accurate vibrational assignments, bond lengths, bond angles, atomic charges and HOMO-LUMO of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole using DFT/B3LYP method. The optimized geometry obtained from DFT calculation was then used to perform NBO analysis. A systematic study on vibrational spectra and structure of the compound is carried out. Numerous reports have been made citing the success of DFT compound to conventional methods in computing molecular and chemical properties such as geometries, vibrational frequencies thermodynamical properties11

 

Experimental Details:

The regents and solvents used were purchased from Sigma Aldrich, Merck specialties Pvt. Ltd and Hi-media laboratories Pvt. Ltd. The compound 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole was prepared according to the following method. A solution of 1-aryl-3-(N, N-dimethyl) thiourea (1mmol) in DMF(2ml) was added to a solution of 2-(2-bromoacetyl)-N-methyl (-6-methylbenzimidazole (0.254g, 1mmol) which was prepared from 2-(1-hydroxyethyl)benzimidazole in DMF (2ml).The reaction mixture was stirred well and triethyl amine (0.15 ml, 1mmol)was added. The reaction mixture was warmed at 80-85°C for 5 minutes. It was then cooled and poured in to ice-cold water with constant stirring. A yellow precipitate thus obtained was filtered, washed with water and dried. The crude product was crystallized from ethanol-water (2:1) to give yellow crystalline solid.

 

Computational Details:

The DFT computation of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methyl benzimidazole has been performed using Gaussian ’09 program package at the Becke-3Lee-Yang-Parr(B3LYP) level with standard 6-31G basis set. The optimized structural parameters are used in the vibrational frequency calculations at DFT level. At the optimized geometry of the title molecule no imaginary frequency modes are obtained, so there is a true minimum potential energy surface is found.

 

The assignments of the normal modes of vibration for the titled compound have been made by visual inspection of the individual mode using the Gauss view software12. All the calculations were done for the optimized structures in gas phase. The methodology involved in density functional theory includes electron correlation and hence is the desired method for obtaining theoretical charge density in molecules. DFT, over the years has become a practical tool for calculating charge-density distributions. Since its adaptability to high- speed computers in easy.

 

RESULTS AND DISCUSSION:

Optimized Geometry:

Optimized geometry was subjected to optimization in the ground state. The optimized structural parameter was calculated by DFT/B3LYP -6-31G basis set is listed in Table-1 in accordance with the atom numbering scheme given in figure-1.


 

 

Figure 1: Optimized geometrical structure of 2-(2-phenylaminothiazol-5-oyl)-1-methyl-6- methylbenzimidazole.

 

Table.1: Optimized geometrical parameters of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole at B3LYP with 6-31G level.

Parameters

Bond lengths(Å)

Parameters

Bond angles(◦)

Parameters

Dihedral angle(◦)

Calculated

Calculated

Calculated

N1-C2

C2-N3

N3-C4

C4-C5

C5-C6

C6-C7

N3-C8

N1-C9

C2-C10

C10-S11

S11-C12

C12-N13

N13-C14

C14-C15

C12-C16

C16-C17

C17-C18

C18-C19

C19-C20

S20-C21

C12-N22

C10-O23

C5-H24

C6-H25

C7-H26

N22-H27

C21-H28

C17-H29

C20-H30

C18-H31

C19-H32

N3-C33

C33-H34

C33-H35

C33-H36

C4-H37

C37-H38

C37-H39

C37-H40

C14-H41

1.3343

1.4000

2.5937

1.3984

1.4147

1.385

1.4004

1.3849

1.475

2.803

1.8384

1.3222

1.377

1.376

2.5129

1.408

1.394

1.400

1.399

1.3986

1.3594

1.2649

1.0860

1.084

1.0835

1.010

1.080

1.087

1.085

1.085

1.084

1.471

1.090

1.090

1.084

1.513

1.097

1.097

1.093

1.079

N1-C2-N3

C2-N3-C4

N3-C4-C5

C4-C5-C6

C5-C6-C7

C2-N3-C8

C2-N1-C9

N1-C2-C10

C2-C10-S11

C10-S11-C12

S11-C12-N13

C12-N13-C14

N13-C14-C15

S11-C12-C16

C12-C16-C17

C16-C17-C18

C17-C18-C19

C18-C19-C20

C19-C20-C21

S11-C12-N22

C2-C10-O23

C4-C5-H24

C5-C6-H25

C6-C7-H26

C12-N22-H27

C20-C21-H28

C16-C17-H29

C19-C20-H30

C17-CH18-H31

C18-C19-H32

C2-N3-C33

N3-C33-H34

N3-C33-H35

N3-C33-H36

N3-C33-H37

C4-C37-H38

C4-C37-H39

C4-C37-H40

N13-C14-H41

112.1955

119.77

138.34

123.63

120.70

106.56

105.98

124.42

155.14

113.69

114.55

112.13

117.50

144.94

141.37

120.22

119.29

121.25

119.09

120.20

117.84

119.04

122.55

115.83

115.83

121.47

119.82

119.91

119.48

120.28

125.79

110.31

110.30

109.01

102.44

112.43

112.43

109.77

120.02

N1-C2-N3-C4

C2-N3-C4-C5

N3-C4-C5-C6

C4-C5-C6-C7

N1-N2-N3-C8

N3-C2-N1-C9

C9-N1-C2-C10

N1-C2-C10-S11

C2-C10-S11-C12

C10-S11-C12-N13

S11-C12-N13-C14

C12-N13-C14-C15

C10-S11-C12-C16

S11-C12-C16-C17

C12-C16-C17-C18

C16-C17-C18-C19

C17-C18-C19-C20

C18-C19-C20-C21

C10-S11-C12-N22

N1-C2-C10-O23

N3-C4-C5-H24

C4-C5-C6-H25

C5-C6-C7-H26

S11-C12-N22-H27

C19-C20-C21-H28

C12-C16-C17-H29

C18-C19-C20-H30

C16-C17-C18-H31

C17-C18-C19-H32

N1-C2-N3-C33

C2-N3-C33-H34

C2-N3-C33-H35

C2-N3-C33-H36

C2-N3-C4-C37

N3-N4-C37-H38

N3-C4-C37-H39

N3-C4-C37-H40

C12-N13-Cl4-H41

-0.0

0

-0.0

0

-0.0

-0.0

179.98

0

-0.0

0

-0.0

-0.0

-179.98

0

179.99

0

0

0

-179.99

-179.97

179.98

-179.99

179.99

-0.0

-179.99

-0.0

-180.0

-180.0

-180.0

179.99

-59.42

59.26

179.9

179.99

-60.77

60.92

-179.92

-179.99

 


The molecule contains Benzimidazole ring, phenyl ring, amino group, methyl group. Geometry optimization can usually locate transition structures. The potential energy surface (PES) specifies hence the energy of a molecular system varies with small changes in the structure. A potential energy surface is a mathematical relationship linking molecular structure and the resultant energy. For a diatomic molecule, it is a two dimensional plot with the inter- nuclear separation Vs potential energy at that bond distance, producing a curve. For larger systems the surface has as many dimensions as there are degrees of freedom within the molecule. Generally a non-linear N atomic molecule has 3N-6 degrees of freedom. The title molecule has 41 atoms and it has 117 degrees of freedom. The optimized bond lengths and bond angles of the title compound are given in Table-1. From the experimental values of literature C-C single bond length is 1.4 A°, C-O length is 1.26 A°, C-N single bond length is 1.40 A°, C-H bond length is 1.08 A°, for benzimidazole ring in title compound and in phenyl ring the C-C bond length is 1. 39 A°, C-H bond length is 1.08 A°, C-N bond length 1.3 A°. The C-C bonds in benzimidazole ring and phenyl ring are not of the same length. Substitution of methyl group leads to some changes of bond angles in the benzimidazole ring. The B3LYP bond lengths are close to the literature data due to slightly exaggerated electron correlation effect.

 

Vibrational Analysis:

In order to obtain the spectroscopic signature of the title compound we performed a frequency calculation analysis13. The scaling factor 0.96 is used for getting theoretical vibrational frequencies. After applying the scaling factor, the theoretically computed wave numbers are in good agreement with literature values. The stimulated IR spectrum of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methyl benzimidazole is shown in figure-2. 


 

 

Figure-2: Stimulated IR spectrum of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methyl benzimidazole.


C-H Vibrations:

The aromatic C-H stretching vibrations are normally found between 3100 to 3000 cm-1 14-16. In this region the bands are not affected appreciably by the nature of the substituents. In benzimidazole ring the C-H stretching mode is observed at 3057cm-1. Our theoretically computed wave number for benzimidazole ring falls within the range of 3092 cm-1 to 3053 cm-1. Thus there is a good agreement between theoretically computed values and literature values. The C-H in plane bending vibrations is observed in the region 1396 cm-1-1320 cm-1. The C-H out of plane bending vibrations of the title compound is well identified within their characteristic region 1027 cm-1-886 cm-1.

 

Ring Vibration:

The C=C vibration are more interesting if the double bonds are in conjugation with the ring. The ring carbon-carbon stretching vibrations occur in the region 1600-1300 cm-1 and are due to the stretching and contraction of all the bands contained in the ring and also the interactions between these stretching vibration include ring breathing modes near 1000v cm-1. In the present work the frequencies are observed at 1601 cm-1. The ring puckering modes usually occur at 645 cm-1- 566 cm-1.

 

Methyl group Vibration:

The C-H methyl group stretching vibrations are generally observed in the range 3580-3095 cm-1 for the assignments of CH3 group frequencies 9 fundamental vibrations can be associated to CH3 group namely CH3 sym stretching, CH3 asymmetric stretching, CH3 in plane scissoring, CH3 out plane scissoring, CH3 in plane bending, CH3 out plane bending, CH3 in plane twisting, CH3 out plane twisting, CH3 torsion modes, for the title compound the methyl in plane bending modes occur in the range 1411 cm-1. The C-H out of plane bending occurs at 837 cm-1 and 764 cm-1. The assignments are in line with the literature values.

 

N-H Vibration:

The hetero aromatic molecule containing an N-H group and its stretching absorption occur in the region 3500-3220 cm-1. Primary amines examined in dilute solution display two weak absorption bands one near 3500 cm-1 and the other near 340 cm-1. These bands represent respectively the asymmetric and symmetric N-H stretching modes17. The theoretically computed wave numbers for the above vibration are observed at 3484 cm-1 represents the asymmetric N-H stretching mode. Primary aromatic amines normally absorb at 1615-1580 cm-1. In the present work, the theoretical calculations indicate the scaled frequency value at 1601 cm-1 is assigned to N- H in plane bending vibration. The presence of aromatic N-H out of plane bending vibration are appeared in the region 767-670 cm-1 for the title compound the theoretical calculation indicates the scaled frequency value at 647 cm-1 is assigned to N-H out of plane bending vibration.

 

Carbonyl group Vibrations:

The carbonyl group is important and its characteristic frequency has been extensively used to study a wide range of compounds. The C=O stretching vibration is observed in the region of 1850-1560 cm.-1 18 . The intensity of these bands can increase due to conjugation or formation of hydrogen bands. The lone pair of electron on oxygen also determines the nature of the carbonyl group. In our present study the theoretically computed wave number for C=O stretching vibration occur at 1556 cm-1.

 

C-S Vibrations:

The C-S stretching vibration is expected in the region 710-690 cm-1 19. In the present work, the frequency is observed at 687 cm-1 which is in line with literature value.

 

C-N Vibration:

The identification of C-N vibrations is a difficult task20, since the mixing of vibrations is possible in this region 1600-1500 cm-1. In the present work the theoretical calculation (benzimidazole ring) indicates the scaled frequency at 1556 cm-1 which is assigned to C-N vibration.


 

IR absorption frequency of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole

Modes

Exp IR frequency

Frequency  (unscaled)

Frequency (scaled)

Intensity

Vibrational Assignments

117

116

115

114

113

112

111

110

109

108

107

106

105

104

103

102

101

100

99

98

97

96

95

94

93

92

91

90

89

88

87

86

85

84

83

82

81

80

79

78

77

76

75

74

73

72

71

70

69

68

67

66

65

64

63

62

61

60

59

58

57

56

55

54

53

52

51

50

49

48

47

46

45

44

43

42

41

40

39

38

37

36

35

34

33

32

31

30

29

28

27

26

25

24

23

22

21

20

19

18

17

16

15

14

13

12

11

10

9

8

7

6

5

4

3

2

1

3447

 

 

 

 

 

 

3167

 

 

 

 

 

 

 

 

 

1604

 

 

 

1580

 

 

1486

 

 

 

 

 

 

1418

 

 

 

 

 

1317

 

 

 

 

 

 

 

 

 

1243

 

1104

 

 

 

 

 

 

 

 

1081

 

 

 

 

 

892

 

 

 

 

 

858

 

 

 

757

 

 

 

 

 

 

 

 

 

 

 

 

 

 

693

 

 

 

 

 

 

3621.79

3297.60

3270.30

3244.22

3238.44

3227.08

3214.52

3210.27

3199.40

3192.74

3178.10

3174.53

3128.93

3107.85

3095.92

3038.34

1665.02

1664.79

1661.45

1635.28

1618.15

1600.59

1583.17

1557.88

1553.61

1550.51

1546.87

1538.86

1536.34

1524.35

1524.18

1506.84

1503.88

1467.37

1447.47

1416.58

1398.39

1396.66

1373.95

1356.72

1329.76

1320.66

1305.88

1300.74

1280.03

1249.47

1239.60

1223.51

1216.92

1171.61

1163.95

1138.07

1133.27

1125.72

1122.0

1097.93

1082.38

1067.64

1039.59

1033.96

1009.20

1004.28

1001.24

978.91

947.67

930.25

921.25

870.61

867.25

843.05

834.45

820.37

804.29

794.51

780.66

776.19

721.40

714.71

672.65

671.31

657.97

653.76

643.95

621.43

614.18

602.63

598.12

588.76

541.38

527.44

507.42

506.75

494.84

428.95

390.54

357.37

342.18

332.01

316.06

303.84

298.39

278.49

274.15

234.26

218.10

205.77

167.30

123.35

121.23

119.27

96.26

67.96

47.14

44.58

37.59

24.49

16.35

3484.0 0

3172.20

3146.02

3120.90

3115.30

3104.40

3092.36

3088.27

3077.82

3071.41

3057.33

3053.89

3010.00

2988.90

2978.20

2922.88

1601.70

1601.50

1598.30

1573.13

1556.60

1539.70

1523.00

1498.68

1494.57

1491.59

1488.0

1480.38

1477.95

1466.42

1466.26

1449.58

1446.73

1411.60

1392.45

1362.74

1345.25

1343.58

1321.73

1305.16

1279.22

1270.47

1256.25

1251.31

1231.38

1201.99

1192.49

1177.01

1170.67

1127.08

1119.28

1094.82

1090.20

1082.94

1079.36

1056.20

1041.24

1027.06

1000.0

994.66

970.850

966.11

963.19

941.71

911.65

894.90

886.24

837.52

834.29

811.01

802.74

789.19

773.72

764.31

750.99

746.69

693.98

687.55

647.08

645.80

632.96

628.91

619.47

597.81

590.84

579.73

575.39

566.38

520.80

507.39

488.13

487.49

476.03

412.64

375.69

343.78

329.17

319.39

304.04

292.29

287.05

267.90

263.73

225.35

209.81

197.95

160.94

118.66

116.62

114.73

92.60

65.37

45.34

42.88

36.16

23.55

15.72

66.9156

9.0810

12.9602

14.7977

15.2929

36.1869

23.5866

28.4396

0.02415

14.2291

13.9492

7.2241

24.4197

16.4869

19.7193

32.3563

16.1779

6.7042

120.8584

9.4612

159.2258

511.3247

229.6870

47.8009

143.6085

80.5184

1.8989

904.3572

253.7118

38.3113

15.2308

363.3581

78.3816

0.3155

25.1985

55.5925

101.9328

145.6871

11.22

4.8989

21.9024

4.2873

43.1703

19.4752

101.9092

552.9616

31.9497

2.9575

4.5103

0.1420

99.3152

8.7425

24.6808

64.5515

3.3098

3.4626

0.7062

3.9370

3.0465

0.9468

10.3459

0.4111

0.3717

6.740

12.8352

0.1108

43.6075

0.5509

5.9584

0.6029

96.4128

30.3619

0.1207

90.3441

10.2044

25.4908

24.2454

3.9633

9.1954

7.9327

67.4585

0.8312

7.7683

8.01511

0.0896

33.3541

0.0896

13.5640

8.7584

20.8107

6.4325

4.9264

0.1951

0.0024

9.5536

3.0948

6.3524

0.0723

3.2719

6.7296

0.7149

0.0576

5.5732

0.2353

2.8753

13.3144

2.3508

0.4003

1.0591

1.2109

1.1046

0.4923

0.4315

0.4853

0.6302

0.2040

2.2726

N22-H27 asy.str.

C14-H41 asy.str.

C21-H28 asy.str.

C33-H36 asy.str.

C7-H26, C6-H25 asy.str.

C18-H31, C28-H30, C19-H32 asy.str.

C6-H25, C7-H26, C5-H24 asy.str.

C18-H31, C17-H29 asy.str.

C20-H30 asy.str.

C5-H24, C6-H25 asy.str.

C5-H24, C6-H25 asy.str.

H34-H35 asy.str.

H38-H39-H40 asy.str.

H34-H35-H36 sym.str.

C18-H31 bend(roc)

H38-H39-H40 sym.str.

Phenyl ring puckering, C-C sym.str.

N22-H27 (bend), C-H bend(wagg.)

C-H bend(wagg.)

C-H bend(wagg.)

C10-O(bend)

N-H bend(roc.)

N-H bend(roc.)

C5-H24 bend(wagg.), C6-H25 bend(wagg.)

Phenyl ring wagging

Methyl group roc.

Methyl group roc.

N22-H27, C19-H32, C20-H30, C18-H31(opb)

Methyl group vib., phenyl ring vib.

Benzimidazole ring vib.

Benzimidazole ring vib.

C5-H24, C6-H25, C7-H26 bend(twist.)

Methyl group vib.

Methyl group sci.

Methyl group twist.

Benzimidazole ring vib., N22-H27 bend(wagg.)

Phenyl ring vib.

C5-H24, C7-H26, C6-H25 bend(wagg.)

C-H in Phenyl ring twist.

Phenyl ring bend(wagg.)

Benzimidazole ring, Phenyl ring vib.

Phenyl ring puckering

C14-H37, N22-H27 bend(roc.)

C5-H24, C6-H25, C7-H26 bend(wagg.)

Phenyl ring vib., C6-H25 ipb, Methyl group vib.

C14-H37 bend (wagg.), Phenyl ring vib.

Phenyl ring vib. Methyl group vib.

Phenyl ring vib. Methyl group vib.

H37-C14bend(twist.)C5-H24, C6-H25 bend(wagg.)

Phenyl ring opb, Methyl group twist.

C19-H32, C18-H31, C20-H30 bend(roc.)

Benzimidazole ring vib., C7-H26 bend(wagg.)

Phenyl ring opb

C15-C14-H37 bend(sciss.), Methyl group opb

C15-C14-H37 bend(sciss.), Methyl group opb

Methyl group bend(roc, )

Phenyl ring puckering

C-H(wagg.) Phenyl ring only

C-H (twist.) Phenyl ring only

Phenyl ring puckering

C-H( roc.) in Phenyl ring

C7-H26, C6-H25, C5-H24(opb)

C17-H29, C18-H31, C19-H32, C20-H30(opb)

C14-H37(opb)

C-H(tiwist.)in phenyl ring

C-H(roc.)in phenyl ring

C-H(wagg.)in phenyl ring

C-H(sciss.)in phenyl ring

Phenyl ring puckering

Benzimidazole ring puckering, methyl group (opb)

Phenyl ring puckering

C20-H30, C19-H32, C18-H31(wagg.)

C-O(twist.), N-C(twist.)

C-H(twist.) in phenyl

C20-H30, C19-H32, C18-H31(wagg.)

Benzimidazole ring puckering

C-H(twist.)in phenyl ring

C-S (Sym.str.), phenyl ring vib.

C15-S11(Sym.str.)

C15-S11(Sym.str.), phenyl ring puckering

C-H(roc.)in phenyl ring

Whole molecule vib.

Whole molecule vib. and phenyl ring puckering

N13-C12-H27(opb)

C-C(twist.), N-C(roc.)

C15-S11(sym.str), benzimidazole ring vib.

C33-H34-H35(roc.), C5-H24, C6-H24(sciss.)

Methyl group(opb), phenyl ring puckering

C15-S11(wagg.), N22-H27(roc.)

N22-H27(twist.), C7-H36, C6-H25(Sciss.)

Benzimidazole ring puckering

N22-H27(roc.)

Phenyl ring(bend)

H29-H31(opb), phenyl ring vib.

Phenyl ring vib., C10-O23(opb)

Phenyl ring vib., C10-O23(opb), methyl group(roc.)

Methyl group (opb)

Benzimidazole ring puckering

C10-O23(roc.)

N22-H27, C2-N1, H24-C5, H25-C6(wagg.)

N22-H27, C2-N1, H24-C5, H25-C6(wagg.)

C37-H38, H39-H40(roc.)

N22-H27(wagg.)

C19-H32, C10-O23(twist)

C37-H38(opb)

C10-O23(roc.)

C17-H29, C18-H31(wagg.)

C10-O23(opb)

Methyl group Vib., phenyl ring puckering

Methyl group Vib., phenyl ring puckering

C33-H34-H35-H36(twist.)

C37-H389(opb)

C37-H389(opb)

C37-H389(opb)

C37-H389(opb), phenyl ring vib.

C33-H34(roc.)

C33-H34(roc.)

 


Mulliken Atomic Charges:

Mulliken atomic charge calculation has an important role in the application of quantum chemical calculation of molecular system because of atomic charges effect dipole moment, molecular polarizability electronic structure and more a lot of properties of molecular systems. The bonding capability of a molecule depends on the chelating atoms. Atomic charge distributions were calculated by determining the electron population of each atom as defined by the basis set. To validate the reliability of our results the mulliken population analysis of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methyl benzimidazole has been calculated using B3LYP/6-31G basis set shown in Table 1. The mulliken atomic charge of sulfur carries positive charge (0.429). Nitrogen has a maximum negative charge value of about -0.436. However all the hydrogen atoms exhibit net positive charge. These magnitudes are changing between 0.119 and 0.339.It is worthy to mention that C2, C4, C8, C9, C10, C12, C14, C16 atoms exhibit positive charge, while C5, C7, C15, C17, C18, C19, C20, C21, C33, C37 atoms exhibit negative charge. The presence of large negative charge on nitrogen and oxygen atom and a net positive charge on hydrogen atom may suggest the formation of intramolecular interaction.

 

Mulliken population analysis of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methyl benzimidazole performed at B3LYP/6-31G set.

Atom

Mulliken atomic charges

Atom

Mulliken atomic charges

N1

C2

N3

C4

C5

C6

C7

C8

C9

C10

S11

C12

N13

C14

C15

C16

C17

C18

C19

C20

-0.473

0.456

-0.747

0.186

-0.189

-0.135

-0.789

0.289

0.085

0.325

0.429

0.172

-0.436

0.517

-0.392

0.313

-0.164

-0.135

-0.111

-0.154

C21

N22

O23

H24

H25

H26

H27

H28

H29

H30

H31

H32

C33

H34

H35

H36

C37

H38

H39

H40

C41

-0.294

-0.758

-0.463

0.119

0.126

0.144

0.339

0.192

0.124

0.132

0.130

0.127

-0.223

0.187

0.187

0.156

-0.482

0.159

0.159

0.151

0.215

 


 

Mulliken charge distribution of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole

 


HOMO-LUMO:

The most important orbital’s in a molecule are the frontier molecular orbital’s, called highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). These orbital’s determine the way the molecule interacts with other species. The lowering of the HOMO-LUMO band gap is essentially a consequence of the large stabilization of the LUMO due to the strong electron- acceptor ability of the electron- acceptor group21. Analysis of wave function indicates that the electron absorption corresponds to the transition from the ground to the first excited state and in explained by one electron excitation from HOMO to LUMO. The frontier orbital gap helps to characterize the chemical reactivity optical polarizability and chemical hardness-softness of a molecule. The HOMO-LUMO energy gap is an important value for stability index22. This reflects the chemical activity of the molecule. The lowering of energy gap explains the eventual charge transfer interactions taking place within the molecule. The molecules having large energy gap are known as hard and molecules having small energy gap are known as soft molecules23. The soft molecules are more polarisable than hard ones because they need small energy to excitation. According to koopman’s theorem24, the HOMO energy is related to the ionization potential and the LUMO energy has been used to estimate the electron affinity. The average value of the HOMO and LUMO energies is related to the electronegativity defined by mulliken

 

ƞ= (ƐHOMOLUMO)/2

 

The calculated energy value of HOMO is -0.24488 and for the LUMO is -0.03399, the value of energy separation ∆E is +0.27887 a.u. Hence from the calculations we conclude that the molecule taken for investigation belongs to soft material.


 

 

HOMO

 

LUMO

 

 

 


Dipole Moment:

DFT predicts dipole moment and higher multipole moments. The dipole moment is the first derivative of the energy with respect to an applied electric field. It is a measure of the asymmetry in the given as a vector in three dimensions. The dipole moments are broken down in to X, Y and Z components. The total dipole moment of the title compound is found to be 3.2685 Debye.

 

CONCLUSION:

The geometry of 2-(2-phenylaminothiazole-5-oyl)-N-methyl-6-methylbenzimidazole was optimized in different levels with DFT-B3LYP method using 6-31G basis set. Using the optimized geometry the vibrational frequencies have been found to agree well with the literature reported values. Atomic charge distributions were calculated by determining the electron population of each atom. The lowering of HOMO and LUMO energy gap explains the charge transfer interactions that take place within the molecule.

 

ACKNOWLEDGEMENT:

T.F. Abbs Fen Reji thanks University Grants Commission, New Delhi for Financial Assistance in the form of Major Research project [F.No.41-229/2012 (SR)]. The authors thank NIIST, Trivandrum and CDRI, Lucknow for spectral and analytical data.

 

REFERENCES:

1.      Mokle, S.S.; Sayeed, M.A.; Kothawar, C; Int.J.Chem.Sci., 2004, 2, 96-100.

2.      Hunter, R.F., J. Chem.Soc., 1925, 2023.

3.      Hunter, R.F., J. Chem.Soc., 1925, 2270.

4.      Joshua, C.P.; Rajasekharan, K.N., Chem.Ind., 1974, 750.

5.      Sareen, V.; Khatri, V., Jain, P.; Sharma, K., Indian.J.Chem., 2006, 45B, 1288.

6.      Sawhney, S.N.; Singh, J., Indian J. Chem., 1970, 8, 882.

7.      Silverstein, M.; Clayton, G.; Morill, B.C., Spectrometric identification of organic compounds, Wiley, New York, 1981.

8.      Balachandrana, V.; Karthick, T.; Perumal, S..Nataraj, A., Spectrochem. Acta A, 2012, 92, 137-147.

9.      Krishnakumar, V.; Prabavathy, N., Spectrochem .Acta A, 2009, 72, 743-747.

10.   Patil, V.S.; Padalkar, V.S.; Tathe, A.B.; Gupta, V.D.; Sekar, N., J. Fluores., 2013, 2015, 1019.

11.   Sawant.A.B.; Nirwan, R.S., Int.J. Chem. Struct., 2013, 47-65.

12.   Mitchison, D.A.; Natur. Med.; 1996, 2, 6.

13.   Pattan, S.R., Suresh, C.H., Pooja, V.D., Indian J.Chem., 2005, 44B, 2404.

14.   Frisch, M.J. et al., Gaussian 09, Revision D.01, Gaussian, Inc., Wallingford CT, 2009.

15.   Obot, I.B.; Obi-Egbedi, N.O.; Umoren, S.A., Int. J. Electro. Chem.Sci., 2009, 4, 863-877.

16.   Dong, S.L.; Cheng, X., Acta Cryst., 2012, 65, 518.

17.   Koopmans, T.A., Physica, 1933, 1, 104.

18.   Kuwae, A.; Machida, K., Spectrokim. Acta A, 1979, 35, 841.

19.   Green, J.H.S.; Harrison, D.J., Spectrochim. Acta A, 1976, 32, 1279.

20.   Sinta, S.P.; Chatterjee, C.L., Spectrosc. Lett., 1976, 9(1976)461.

21.   Patel, V.H.; Gandhi, S, A., Indian J.Pure & Appl. Phys., 2011, 49, 227.

22.   John M. Chalmers, Pete Griffiths(Eds), Handbook of Vib. Spectrosc., Wiley, Newyork, 2002, 5.

23.   Rastogi V.K.; Arora, C.B.; Singhal, S.K.; Singh, D.N. Yadav, R.A., Spectrochim. Acta A, 1997, 53, 2505.

24.   Chandra, R.; Singh, A.; Singh, T.P., Asian J. Phys., 1993, 2, 50.

 

 

 

 

Received on 18.09.2018         Modified on 18.10.2018

Accepted on 15.11.2018         © AJRC All right reserved

Asian J. Research Chem. 2018; 11(6): 848-856.

DOI: 10.5958/0974-4150.2018.00149.9